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An evaluation of the 1965.0 epoch, H. O. 1706 World Magnetic 
Variation Chart is performed by direct comparison of charted 
and observed values. Results of this study provide an estimate 
of chart reliability based on the correlation of the computed 
spherical harmonic surface and individual field data values. 
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FOREWORD 


An important function performed by the U. S. Naval Oceanographic 
Office is the periodic publication of the H. O. 1706 World Magnetic Varia- 
tion Chart. The recent edition of this chart was, for the first time, compiled 
from a mathematical model of the geomagnetic field. This report presents the 
results of an evaluation of the chart's reliability and provides insight into the 


challenging task its construction poses. 
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INTRODUCTION 
Background 


The Naval Oceanographic Office in conjunction with the Coast and 
Geodetic Survey, publishes the H. O. 1706 World Magnetic Variation Chart 
every five years. This chart is a smoothed representation of the earth's mag- 
netic field depicted by contours of equal values called isogonic lines. Such 
charts are indispensable, because from them are derived data necessary to 
accurately navigate both military and commercial ships and aircraft. 


For a number of years the cartographic technique used in compiling 
the world chart consisted of manually plotting and contouring smoothed survey 
data extrapolated to epoch. With the publication of the 1965 world chart, 
however, a more sophisticated approach in constructing the chart was 
employed. The Coast and Geodetic Survey, for the first time, constructed 
the world chart from a mathematical model of the geomagnetic field. An 
important step in the derivation was the analytic smoothing of survey data. 
The smoothing produced a mathematical model expressed by a set of spherical 
harmonic coefficients to degree and order 12 (Hurwitz et al., 1966). Chart 
isolines were computed from the mathematical model and automatically plotted 
to produce the world chart. 


Purpose of Report 


Because of the important role magnetic variation plays in navigation, 
its accurate portrayal is a significant goal. To this end, evaluations of 
previous world charts have been made with subsequent improvements incor- 
porated into the charts. A detailed evaluation, however, has not been 
practical prior to the 1965 chart. This has been due, in large part, to the 
formidable task of abstracting representative quantities of data from the 
charts for comparison with observed data. With coefficients of a mathema- 
tical model available for the first time, the difficulty of synthesis can be 
resolved by electronic computers. Utilizing this convenience, a direct com- 
parison of charted and observed values was made and an estimate of the 1965 
World Magnetic Variation Chart's reliability achieved. 


PRELIMINARY PROCESSING 
Division of World Chart 


To perform the evaluation in as much detail as practical, the World 
Aeronautical Chart quad-system was chosen to divide the world into small 


areas. This system produces a total of 1,851 charts, illustrated as Figure 1. 
The numbering of these charts begins in the northern polar area and progresses 
to the southern polar area. For convenience, the area above 84° north was 
considered as one chart, a slight deviation from the conventional system. 
Otherwise, each chart covers 4° of latitude with the range of longitudes 
depending upon the latitude. The projection is Mercator between latitudes 
0° and 4°, Lambert conformal between latitudes 4° and 72°, and polar 
stereographic between latitudes 72° and 90°. Each chart containing obser- 
vational data was studied independently, with the results based on those 

data falling within the chart boundaries. 


Data 


The observed variation data used in this study are identical to the 
data used in compiling the recent world chart. These data represent approxi - 
mately 165,000 values measured at land, sea, and airborne locations since 
1900. Airborne values were reduced to sea level and all data extrapolated 
to epoch by applying the appropriate secular change rate. In addition to the 
observed data, the corresponding charted values were computed from the 
mathematical model. Finally, a residual value or difference between each 
charted and observed data point was calculated. The above information was 
then sorted by geographic position into the appropriate World Aeronautical 
Chart. 


ANALY SIS 


Following the preliminary preparation of the data, a FORTRAN 
program was written to perform the analysis. An IBM 7074 computer executed 
all computations and prepared the results in a format suitable for plotting on 
the CalComp Automatic Plotter. The resulting graphics and associated statis- 
tics are presented in Table Il. Definitions and formulae pertinent to the 
analysis are discussed below. 


Histograms 


A histogram was used to depict the frequency distribution of the 
residual values of variation within each World Aeronautical Chart. The 
horizontal axis was divided into thirty-one intervals with the first and last 
interval regarded as "open-ended," i.e., without specific limits. All other 
intervals have fixed boundaries with a range of 0.5°. A mid-point value 
for each interval is indicated and labeled in 0.5° increments with zero 
representing that point where charted and observed data are in perfect 
agreement. Treating the charts independently, each residual value of 
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variation was sorted with respect to both sign and magnitude into the appro- 
priate interval and a count of residuals falling in each interval tabulated. 
The vertical axis of the histogram is the percentage of the total number of 
residuals contained in each interval and is graduated in steps of 10%. 


Statistics 
1. World Aeronautical Chart 
The histogram and associated statistics for each World Aero- 
nautical Chart (WAC) are presented in the numeric sequence 
shown in Figure 1. WAC numbers are listed consecutively 
except for those gaps created by charts containing no data. 


2. Number of Points 


The total number of data points (NO PTS) within each 
chart is likewise tabulated and presented with each illustra- 


tion, 
3. Mean 
The mean for each chart was then computed and is defined 
as follows: 
Mean (X;) = 2:%e = Xo) : 20%) 
Where 


Xp = mean of the residual variation values 
X_ = charted variation value (H.O. 1706, epoch 1965.0) 


Xo = _ observed variation value (all available observations 
since 1900 extrapolated to epoch 1965.0) 


Xr = residual variation value 


N total number of data points within the chart 


This value represents a type of average that divides a distribution 
into two equal parts in terms of the total values of the residuals. 


4. Root Mean Square 


Using previous definitions, the root mean square (RMS) was 
derived by: 


RMS = Recta fa fe EWEN 1 


This value provides an indication of the degree of scattering or 
dispersion of the residuals about their mean. 


5. Median 

The median of the residuals was computed next and is an 
average that divides a distribution into two equal parts in terms 
| , of the 


of frequencies. It is the value that leaves 50%, ( 5 


number of residuals above and below it. The formula used for 
obtaining this value is: 


W 
Median = L +(*) | 
Where 
: N 
L = upper boundary value of interval where >) occurs or 
is as close as possible without exceeding it 
; a IN 
W = number of values needed from next interval to obtain > 
f = total number of values in next interval 


range of magnitude of the interval (in this case, 0.5°) 
6. 80 Percent Confidence Interval 


Finally, an 80% confidence interval (80 PCT CONF) for the 


median was computed by the following formula: 


80 PCT CONF = Median +1.28 (Standard Error of the Median) 


1.28 = number of standard errors necessary to reach an 80% 
level of confidence 


STANDARD ERROR = 1.2533 (ANE) 


AN 


By establishing a confidence interval on the basis of the sample, one 
can express with 80% confidence the range within which the true 
median lies. 


DISCUSSION 


Obviously, interpretation of the statistics generated for each World 
Aeronautical Chart depends on the size of the sample. As shown in Table I, 


a large percentage of the charts contain insufficient data to permit detailed 
conclusions. 


TABLE | 
DATA DENSITY AND DISTRIBUTION 


wore NORTHERN HEMISPHERE || SOUTHERN HEMISETIERE 
Seuss) vee |e eee ees ee 


Of the 931 charts comprising the southern hemisphere, 814 or over 87% contain 
less than 50 data points. Although the northern hemisphere has been sampled 
somewhat better, 433 or over 47% of the 920 charts still contain less than 50 
data points. Because of the limited number of observations both the charted 
isogonic lines and the statistics of this evaluation must be weighed carefully in 


these poorly sampled areas. Table | illustrates quite vividly the poor data 
distribution and density throughout the world and, therefore, the difficulty 
a cartographer faces in accurately depicting the geomagnetic field. 


In general, the difference between charted and observed data is a 
fair guide to the reliability of the world chart. Accordingly, the mean and 
median for each World Aeronautical Chart was computed fo investigate these 
differences. In each case both magnitude and direction are indicated. The 
sign associated with each value is positive whenever the charted data exceeds 
the observed data. While both values are indices of central tendency, a rela- 
tive advantage of median over mean does exist, i.e., the mean is greatly 
influenced by extreme values, while the median is not. For this reason and 
because the large volume of data used was not edited, the median appears to 
be the more reliable value. The importance of both mean and median in this 
analysis is that they provide a measure of how closely the charted isogonic 
lines and observed data agree. 


In considering the median, it seemed desirable to establish a level of 
confidence. Therefore, an 80% confidence interval based on the computed 
median was derived. This interval represents the range within which the true 
median would lie 80 out of 100 times. A broad interval then, would indicate 
that a "good" estimate is not possible. In this analysis, the median and 80% 
confidence interval are not valid whenever 50% or more of the data lie in the 
"open-ended" intervals of the histogram. This consistency occurs in charts 
1773, 1774, 1815, 1818, 1821, and appears to result from lack of data and/or 
being located in high latitudes. 


The agreement of charted and observed data in the northern hemisphere 
is somewhat better than that of the southern hemisphere. Figure 2 illustrates 
this, depicting those areas where the median exceeds 0.5°. For the northern 
hemisphere, the median in 65 or about 7% of the 904 charts containing data 
exceeds 0.5°. In the southern hemisphere, the median exceeds 0.5° in 221 or 
about 27% of the 824 charts containing data. When considered overall, only 
286 or about 17% of the 1,728 charts containing data have medians exceeding 
0.5°. In light of the sparseness of observed data and the limited detail the 
world chart is designed to show, the mathematical model of the Coast and 
Geodetic Survey would then appear to fit the observed data remarkably well . 


The root mean square, RMS, for each chart was computed to indicate 
the dispersion or scattering of the residuals about their mean. With a large 
sample and a distribution curve which approaches normality, 68% of the resi- 
duals will fall within an interval given by X + RMS, and 95% within 
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X + 2(RMS). While this information is helpful in detecting anomalous areas, 
the differentation of true anomalies from other variables such as diurnal 
variation would be difficult. 


SUMMARY 


An evaluation of the 1965, H.O. 1706 World Magnetic Variation 
Chart was performed by direct comparison of charted values and survey data. 
Utilizing the World Aeronautical Chart quad-system, a mean and median 
value was computed for each chart providing an accuracy estimate of the 
isogonic lines appearing on these charts. Based on these values it was 
observed that the mathematical model used in constructing the 1965 world 
chart generally agreed with observed data to within 0.5°. An 80% confi- 
dence interval and root mean square were also compiled for each chart and 
provide indices of reliability and dispersion respectively. 


CONCLUSIONS 


It is generally conceded that present world charts are probably accur- 
ate to about 0.5° for the United States and Europe, but more in error in 
other areas, particularly at high latitudes (Chernosky et al., 1965). From 
the results of this study we have shown that 83% of the World Aeronautical 
Charts investigated have a median of less than 0.5°. This indicates that 
previous estimates of chart accuracies, while correct, are perhaps conserva- 
tive. We must realize, however, that a great deal of refinement to the world 
chart remains. 


Improvement of future charts depends, in large part, on the acquisition 
of additional data. As this study shows, over 67% of the areas considered are 
represented by less than 50 data points. This is particularly true for most of 
the southern hemisphere and the polar areas. Although many of these areas 
are inaccessible at present, additional data would result in considerable 
improvement of the world chart. 


The mathematical model derived by the Coast and Geodetic Survey 
is, for all practical purposes, as refined as is warranted by the available 
observed data. It is concluded, therefore, that the recent world chart 
drawn exactly from this model is equally as realistic. 
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